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Abstract. We have used Spitzer and its IRAC camera to search for the transit of the super-Earth HD 40307b. The transit- 
ing nature of the planet could not be firmly discarded from our first photometric monitoring of a transit window because of 
the uncertainty coming from the modeling of the photometric baseline. To obtain a firm result, two more transit windows 
were observed and a global Bayesian analysis of the three IRAC time series and the HARPS radial velocities was performed. 
Unfortunately, any transit of the planet during the observed phase window is firmly discarded, while the probability that the 
planet transits but that the eclipse was missed by our observations is nearly negligible (0.26%). 
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1. Introduction 

Transiting extrasolar planets are and will remain key objects 
for our study and understanding of the large planetary popula- 
tion hosted by our galaxy. Except for the planets of our own 
solar system, transiting exoplanets are the only ones with accu- 
rate estimates of mass, radius, and, by inference, constraints 
i on internal composition. Furthermore, their special geomet- 
rical configuration gives us the opportunity to study directly 
their atmospheres without the challenging need to spatially re- 
solve their light from that of their host star. Occultation pho- 
tometry and spectroscopy has been used to map the thermal 
emission of more than a dozen highly irradiated giant plan- 
ets and to assess their atmospheric composition, thermal gradi- 
ent and albedo (see, e.g., Deming 2009). More demanding in 
terms of signal-to-noise ratio (SNR), transit transmission spec- 
troscopy has been used to detect atomic and molecular features 
in the atmospheres of the two giant planets HD 209458b and 
HD 189733b (e.g., Charbonneau et al. 2002, Vidal-Madjar et 
al. 2003, Swain et al. 2009). These transit and occultation mea- 
surements have opened a new field of astronomy: exoplanetol- 

ogy- 
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One of the next major steps in this field will be the first 
atmospheric characterization of a terrestrial extrasolar planet. 
Radial velocity (RV) surveys have opened the road to this goal 
by revealing a population of planets of a few Earth masses 
(e.g., Udry et al. 2007, Mayor et al. 2009a, 2009b), the so- 
called 'super-Earths'. Last year, the CoRoT space mission de- 
tected the first of these super-Earths that transits its parent star: 
CoRoT-7b (Leger et al. 2009, Queloz et al. 2009). With a mass 
M p = 4.8 + 0.8 M e and a radius R p = 1.7 + 0.1 R e , CoRoT- 
7b has an average density of 5.6 + 1.3 g cirT 3 , similar to the 
Earth's. Still, its structure is quite uncertain: its measured mass 
and radius are consistent with a rocky planet, but also with 
a water-rich planet that could contain up to 40% of ice wa- 
ter (Valencia et al. 2009). Due to the relative faintness of the 
host star (V=11.7, K=9.8) and the small planet-star size ratio 
(3.4 x 10~ 4 ), a thorough characterization of this planet is unfor- 
tunately not possible with existing or planned instruments. 

The detection of a second transiting super-Earth was re- 
cently announced by the MEarth Project (Charbonneau et al. 
2009). This new planet orbits around the nearby M-dwarf 
GJ 1214 (M4.5V, d = 13 pc). Its mass is M p =6.5 + 1 M @ , 
similar to the one of CoRoT-7b. Nevertheless, GJ 1214b has a 
much larger radius than CoRoT-7b: R p = 2.68 + 0.13 R @ . The 
resulting average density of 1 .9 + 0.4 g crrr 3 suggests the pres- 
ence of a significant gas component. GJ 1214b is thus more 
a kind of 'mini-Neptune' than a massive rocky planet, but its 
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actual composition remains uncertain, as shown by Rogers & 
Seager (2009). GJ 1214b is a favorable target for atmospheric 
measurements able to provide valuable insights into its com- 
position and history. In particular, the future JWST telescope 
(Gardner et al. 2006) will have the potential to thoroughly char- 
acterize the atmosphere of this super-Earth, mainly because of 
the infrared brightness (K = 8.8) and small size (0.21 R Q ) of its 
host star. 

The future studies made possible by the detection of 
GJ 1214b outline how important are the few low-mass plan- 
ets that must transit nearby bright stars. Due to their need to 
observe thousands of stars at once, most transit surveys do 
not target the brightest and most nearby stars. On the con- 
trary, Doppler surveys target a large number of very bright 
stars. In this context, it is not surprising that the two best stud- 
ied hot Jupiters, HD 209458b and HD 189733b, were first de- 
tected by RV measurements and caught afterwards in transit 
(Charbonneau et al. 2000, Henry et al. 2000, Bouchy et al. 
2005). This was also the case for the hot Neptune GJ 436b 
(Butler et al. 2004, Gillon et al. 2007b), an excellent target too 
for detailed follow-up atmospheric studies (see, e.g., Deming 
et al. 2007, Demory et al. 2007). Searching for the transits of 
the low-mass planets detected by RV measurements is thus an 
obvious method to detect transiting super-Earths suitable for a 
thorough atmospheric characterization. Doppler surveys have 
now detected enough nearby low-mass planets to make highly 
probable that at least some of them transit their parent stars. 
In particular, our HARPS Doppler survey (Mayor et al. 2003) 
has now detected more than 40 low-mass planets (Lovis et al. 
2009). Among them are several announced hot Neptunes (see, 
e.g., Santos et al. 2004, Bouchy et al. 2009) and super-Earths 
(see, e.g., Udry et al. 2007, Forveille et al. 2008), but also many 
more firmly confirmed planets that are waiting for a final anal- 
ysis or some more RV measurements to fully characterize their 
system. 

Differential photometric precisions of 0. 1 % are now rou- 
tinely achieved from the ground, making rather easy the de- 
tection of the transit of a hot Jupiter. Due to the small size of 
its host star, the transit of the hot Neptune GJ 436b could also 
be detected from the ground (Gillon et al. 2007b). For most 
of the low-mass planets detected by HARPS around solar-type 
stars, a transit detection requires a much more challenging pho- 
tometric precision, ~50 ppm or better on a time scale of a few 
hours (corresponding to the mean transit duration for a close-in 
planet), and for very bright stars. From the ground, such a pre- 
cision is presently out of reach. While the photon noise would 
not be a problem even for relatively small telescopes, any atmo- 
spheric instability creates a correlated noise with a magnitude 
larger than the one cited above. Differential photometry helps 
reducing the amplitude of this correlated noise, and ground- 
based sub-mmag transit light curves nearly free of correlated 
noise have been obtained for a few transiting planets (see, e.g., 
Gillon et al. 2009, Johnson et al. 2009, Winn et al. 2009), but 
for rather faint stars (V > 10) with several comparison stars 
of similar brightness nearby (at a few arcmin at most), mak- 
ing possible an excellent correction of atmospheric effects. For 
bright nearby stars, such an optimal configuration is not avail- 
able. Furthermore, the transit ephemeris of the very low-mass 



planets detected by HARPS have generally a 2-cr probability 
interval of a dozen of hours or more. This is not only because of 
the small amplitude of the RV signal due to the planet, but also 
to the RV low-frequency noise of the host star and the appar- 
ently large occurence of multiple low-mass planetary systems 
(Lovis et al. 2009). Due to the day-night cycle and to the fact 
that only a fraction of some nights is suitable for precise pho- 
tometry of the target, ground-based telescopes are definitely not 
the best option for such a transit search, and a space-based in- 
strument is required. 

The needed space telescope has to be able to perform 
exquisitely precise photometry for bright stars but also to moni- 
tor continuously the same star during dozens of hours. We have 
concluded that the best choice for such a program would be the 
Spitzer Space Telescope (Werner et al. 2004). Due to its helio- 
centric orbit, it can monitor most of the stars for several weeks 
during their visibility windows, but it has also demonstrated 
at many instances its excellent photometric potential, detecting 
signals with an amplitude of few hundreds of ppm (see,e.g., 
Beaulieu et al. 2007). We have thus set up a Spitzer program de- 
voted to the search of the transits of HARPS low-mass planets. 
This program is presently divided in two Spitzer sub-programs: 
a cycle 5 DDT program (ID 495) devoted to only the planet 
HD 40307b (Mayor et al. 2009a, hereafter M09), and a 100 
hours cycle 6 DDT program (ID 60027) devoted to ~ 10 other 
low-mass planets. 

This paper presents our results for HD 40307b. This super- 
Earth (M p sin i = 4.2M ffi ) orbits around a bright (V = 7.2, K = 
4.8) nearby (12.8 pc) K2-dwarf. Because of its proximity to 
its host star (a = 0.047 AU), its geometric transit probability is 
~ 7%, making it a good candidate for a transit search. Section 2 
presents the data used in this work, including the Spitzer obser- 
vations and their reduction. Our data analysis is presented in 
Sect. 3, and our results are presented and discussed in Sect. 4. 
We give our conclusions in Sect. 5. 

2. Data 

2.1. Radial velocities 

We used the 135 HARPS measurements presented in Mayor 
et al. (2009a) to estimate the HD 40307b transit windows for 
our Spitzer observations. We also used these data in our final 
global analysis aiming to assess the transit probability of the 
planet (see Sect. 3). 

2.2. Spitzer IRAC photometry 

HD 40307b was first observed with the Infra-Red Array 
Camera (IRAC, Fazio et al. 2004) of Spitzer on 31 October 
2008 from 01h02 UT to 09h40 UT, corresponding to a ~ 2- 
cr transit window of the planet as predicted by our RVs. With 
magnitude K = 4.8 and K2V spectral type, HD 40307 is a 
very high signal-to-noise ratio (SNR) target for all IRAC bands. 
Considering the substantial experience of our team with Spitzer 
photometry at 8 fim (e.g., Deming et al. 2007, Gillon et al. 
2007a, Demory et al. 2007), we conservatively decided to ob- 
serve in this channel (SiAs detector) and to use the established 
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technique of continuous staring in non-dithered subarray mode. 
We choose to use an exposure time of 0.32s, the largest one 
for which the star would not be saturated on the detector. The 
Spitzer flux density estimator toofl predicts a flux density of 
806 mJy for HD 40307b at 8 pm, which translates into a SNR 
of ~5100 for a 1 minute integration when taking account the 
photon, read-out and background noises, and the instrumental 
throughput correction^. 

Values for the stellar and planet sizes are needed to estimate 
the SNR that should be expected for the detection of a putative 
transit. Using 7g ff = 4977 + 59 K, logg = 4.47 + 0.16 and 
[Fe/H] = -0.31 + 0.03 (M09), the theoretical stellar calibra- 
tion obtained by Torres et al. (2009) from well-constrained de- 
tached binary systems predicts a stellar mass of 0.78+0.06 M Q , 
in good agreement with the estimation of Sousa et al. (2008), 
0.77 + 0.05 Mq. The parallax of this nearby star was precisely 
determined by Hipparcos (ESA 1997) to n = 77.95 + 0.53 mas. 
From this value, the Hipparcos V magnitude (V - 7.17, ESA 
1997), the 2MASS K magnitude (K = 4.79, Skrutskie et al. 
2006), the spectroscopic parameters from M09 and the bolo- 
metric correction calibration presented by Masana et al. (2006), 
we deduce a radius of 0.68 + 0.02 R Q for the star (correspond- 
ing to a luminosity L„ = 0.26 + 0.01 L ). Depending on the un- 
known planetary composition, the radius of HD 40307b could 
lie between 1 . 1 up to more than 5 R @ (Seager et al. 2007), the 
smallest value corresponding to an unlikely pure iron composi- 
tion. Assuming R„ = 0.68+0.02 R Q , R p - 1 . 1 R® translates into 
a transit depth of 220 ±13 ppm and a maximum transit duration 
of ~ 135 min. Assuming an actual duration of 30 min (and thus 
a large impact parameter), the resulting expected SNR on the 
transit detection is ~ 6. Thus we should have been able to de- 
tect the transit of the planet for any plausible composition and 
for most of the possible impact parameters. 

The data of the run are constituted of 1297 sets of 64 
individual subarray images. Two of us (DD and MG) per- 
formed an independent preliminary reduction of these images. 
Aperture photometry was used by both of us on the Basic 
Calibrated Data (BCB) provided by the standard Spitzer re- 
duction pipeline. After a rejection of the discrepant fluxes, the 
measurements were averaged for each set of 64 images, lead- 
ing thus to 1297 measurements in the final light curves. The 
systematic effect known as the 'ramp' (see, e.g., Knutson et 
al. 2007) is extremely sharp at the beginning of the run (see 
Fig. 3). In his preliminary analysis, DD simply cut off the first 
three hours and fitted a straight line to the remainder. Figure 1 
shows the residuals of this fit vs HJD. A tiny drop in bright- 
ness is visible from -2454770. 70 HJD. Some tests were done 
to assess its significance. A transit curve having a fixed du- 
ration - calculated for a central transit - was constructed and 
fitted as a function of center time and transit depth. The best 
fit is overplotted in Fig. 1. From the precision on its depth and 
on the out-of-transit level, its significance is ~ 3cr. It is thus a 
marginal detection. 50,000 bootstrap trials were done wherein 
data were permuted and the maximum amplitude of each set 
was found, varying the center time. Only 1 -percent of the trials 

1 http://ssc.spitzer.caltech.edu/tools/starpet/ 

2 http://ssc.spitzer.caltech.edu/documents/som/ 
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Fig. 1. HD 40307 Spitzer 8 pm photometry (first run) obtained 
after reduction by DD. The first 3 hours have been cut off and a 
line was fitted in the photometry to correct for the 'ramp'. The 
best-fitting central transit model is overimposed. See text for 
details. 

produced an amplitude larger than the one of the structure de- 
tected. Varying the aperture, this fraction ranged from zero to 
2%, so a 98% confidence could be assigned to the transit-like 
structure under the simple line model for the baseline. 

MG cut off the first hour of data and obtained a light curve 
that shows a similar drop in brightness (see Fig. 2). Its signif- 
icance was tested with a Markov Chain Monte Carlo method 
(see Sect. 3). Two Markov chains were performed. For each 
of them, the free parameters were the transit depth, the impact 
parameter, the central timing and the parameters of the base- 
line model (see Sect. 3 for details). The stellar mass and radius 
were fixed to M, = 0.78 + 0.06 M and R t = 0.68 + 0.02 
R G . For the first chain, the initial transit depth was set to 300 
ppm, while it was set to zero for the second chain. The re- 
sults for both chains were analyzed together to obtain the pos- 
terior distribution for all fitted parameters. The solutions for 
both chains are fully compatible, the deduced transit depth be- 
ing dF = 170^2 ppm. The significance based on this analy- 
sis is thus comparable to the one deduced from the analysis 
performed by DD. The deduced transit depth translates into a 
planet radius R p = 0.97+Q2 3 /?®. Such a small radius could be 
estimated a priori unlikely, because it corresponds to a planet 
denser than expected for a pure iron composition. Still, any un- 
known systematic error in the measurements used to estimate 
the stellar radius could lead to a larger planet. For instance, 
the stellar calibration of Torres et al. (2009) leads to a larger 
stellar radius, 0.86 + 0.17 R . We thus chose not to discard our 
marginal transit detection on the basis that it was shallower than 
expected, and we tried to confirm it with more observations. 

We requested the observation of two more transit windows 
with Spitzer. HD 40307 was again observed with IRAC on 02 
February 2009 from 21h26 to 05h52 UT and on 13 March 2009 
from 16h23 to 00h48 UT. For both runs, 1108 set of 64 indi- 
vidual subarray images were acquired at 8 pm, again with an 
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Fig. 2. HD 40307 Spitzer 8 pm photometry obtained after re- 
duction by MG. The first hour has been cut off and the points 
have been binned per 20 minutes. The best-fitting ramp+transit 
model is overimposed. See text for details. 

exposure time of 0.32s. At this stage, it was decided to obtain a 
final light curve for each run and to perform a global Bayesian 
analysis of these light curves in addition to the HARPS mea- 
surements to assess the probability that HD 40307b was transit- 
ing (see Sect. 3). For the three Spitzer runs, we converted fluxes 
from the Spitzer units of specific intensity (MJy/sr) to photon 
counts, and aperture photometry was obtained for HD 40307 
in each image using the IRAF/DAOPHOT0 software (Stetson, 
1987). A circular aperture with a radius of 6 pixels was cen- 
tered in each image by fitting a Gaussian profile on the target. 
A mean sky background was measured in an annulus extend- 
ing from 16 to 24 pixels from the center of the aperture, and 
subtracted to the measured flux for each image. For each set 
of 64-images, a 3-cr median clipping was used to reject out- 
liers, then the remaining values were averaged and the error 
on the mean was considered as the error on the resulting mea- 
surement. The mean number of measurements rejected per set 
for the three runs together was 1.5. At this stage, a 4-cr clip- 
ping median filtering was used to reject outlier sets. For the first 
run, two measurements were rejected out of 1297. For both the 
second and third run, three measurements were rejected out of 
1 108. The resulting light curves binned with an interval of five 
minutes are shown in Fig. 3 

3. Data analysis 

The three Spitzer 8 pm time series were first binned with an 
interval of five minutes to speed up the analysis. These binned 
photometric time series and the HARPS measurements were 
used as input data in a global analysis aiming to determine 
the probability that HD 40307b transits and, if so, its tran- 
sit and physical parameters. The analysis was performed with 
the adaptative Markov Chain Monte Carlo (MCMC) algorithm 
presented in Gillon et al. (2010). The assumed models were 
based on a star and three planets on a Keplerian orbit about 
their common center of mass (the planet-planet interactions 
were first assessed and revealed to be negligible). The planet 

3 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



HD 40307b had its orbital inclination free and was thus al- 
lowed to transit the star. To model the eclipse photometry, we 
used the photometric eclipse model of Mandel & Agol (2002) 
multiplied by a systematic effect model. For the RVs, a classi- 
cal Keplerian model was used in addition to a linear drift (see 
M09) and a Rossiter-McLaughlin (RM) effect model (Gimenez 
2006) for the RVs obtained during a transit. 

Considering the very red bandpass of our Spitzer photom- 
etry, we assumed no limb-darkening for the photometry. For 
the RVs, a quadratic limb-darkening was assumed. The values 
Mi = 0.613 and m = 0.161 were deduced from Clarets tables 
(2000; 2004) for the stellar parameters presented in M09 and 
for the V-filter, corresponding to the maximum of transmission 
of HARPS. These values were kept fixed in the MCMC. 

The 8 pm IRAC photometry is known to be affected by a 
systematic effect causing the gain to increase asymptotically 
over time for every pixel, with an amplitude depending on 
their illumination history (see e.g. Knutson et al. 2008 and 
references therein). Following Charbonneau et al. (2008), this 
'ramp' was modelled in each IRAC time series as a quadratic 
function of ln(dt): 

A(dt) — a\ + ailn{di) + a^(ln(dt)) 2 + a\dt, (1) 

where dt is the elapsed time since 15 min before the start of 
the run. The linear term a\dt was added to take into account 
that HD 40307 is a slowly rotating K-dwarf star that could be 
variable on a time scale of several weeks, and that such a low- 
frequency variability could translate into a slight slope in the 
Spitzer time series. 

The photometric correlated noise and the RV jitter noise 
were taken into account as in Gillon et al. (2010). For the three 
planets orbiting HD 40307, the jump parameters in the analysis 
were the two Lagrangian parameters e cos a> and e sin u where 
e is the orbital eccentricity and u> is the argument of perias- 
tron, the orbital period P, the K2 parameter characterizing the 
amplitude of the orbital RV signal (see Gillon et al. 2010) and 
the time To for which the true anomaly M = n/2 - a>. For a 
transiting planet, Tq represents the time of minimum light. For 
HD 40307b, the planet/star area ratio (R p /R s ) 2 and the impact 
parameter b' — a cos i/R, were also jump parameters. 

Because of the expected small planet-star area ratio and 
the small rotational velocity of the star (Vsin/ < 1 km.s -1 , 
M09), the amplitude of the RM anomaly in case of transit of 
HD 40307b would be smaller than the typical precision of the 
HARPS measurements, so we kept fixed the projected rota- 
tional spin velocity of the star Vsin / to 1 km.s 1 and the pro- 
jected spin-orbit angle /3 of HD 40307b to deg. 

A value for the stellar mass and radius were drawn at each 
step of the Markov chains from the distributions N(0.18, 0.06 2 ) 
M and A^(0.68, 0.02 2 ) R Q . A uniform prior distribution was 
assumed for each jump parameter. For each Spitzer time se- 
ries, the four parameters of the baseline model were determined 
at each step of the MCMC by least-square minimization (see 
Gillon et al. 2010 for details). It was also the case for the sys- 
temic RV of the star. The planet radius was forced to range 
from 1 to 6 R @ to avoid unrealistic solutions. 

The analysis was divided into two steps. A first Markov 
chain was performed to assess the level of correlated noise in 
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Fig. 3. HD 40307 Spitzer 8 pm photometry binned with an in- 
terval of five minutes. The best-fitting model is overimposed. 
See text for details. 

the photometry and the jitter noise in the RVs and to update the 
measurement error bars accordindly. A value of 0.85 m.s -1 was 
deduced for the RV jitter noise, while /3 re d = [1, 1, 1.14] were 
deduced for the photometric time series (see Gillon et al. 2010 
for details). Five new chains (10 6 steps each) were then per- 
formed using the updated measurement error bars. The good 
convergence and mixing of these five chains was checked suc- 
cesfully using the Gelman and Rubin (1992) statistic. 

4. Results and discussion 

Table 1 shows the median values and 68.3% intervals of the 
posterior distributions for the parameters of the three planets. 
These values are in good agreement with the ones reported by 
M09. Figure 3 shows the three IRAC time series with the best- 
fitting model overimposed. No transit is present for the best- 
fitting solution. The best-fitting RV models are showed with 
the RV measurements in Fig. 4. 

Our global Bayesian analysis leads to a posterior transit 
probability for HD 40307b of 0.3%. If we consider only full 
transits, the posterior probability is only 0.19%. The marginal 
transit detection that we obtained from the first IRAC run is 



S o 




dT [days 



Fig. 4. Period-folded HARPS measurements and Keplerian 
curve for each of the planets, after correction of the effect of 
the two other planets and the drift. Top: HD 40307b, middle: 
HD 40307c, bottom: HD 40307d. 



thus unfortunately not confirmed, and the transiting nature of 
HD 40307b is rejected at ~ 3<x. As shown in Fig. 5, the frac- 
tion of the simulations for which HD 40307b transits concern 
mostly transits that happen outside the observed window. The 
resulting posterior probability that a transit was present in the 
observation window is 0.04%, while no simulation resulted in 
a full transit present in the IRAC runs. The presence of a transit 
in the Spitzer data is thus rejected with a high level of confi- 
dence, while the probability that the planet transits but that the 
eclipse was missed is nearly negligible (0.26%). 

To investigate the effect of the stellar size on this result, we 
performed a second MCMC analysis assuming R t - 0.86+0.17 
R G instead of 0.68 + 0.02 R G (see discussion in Sect. 2). The 
obtained transit probability was 0.63 %, leading to the conclu- 
sion that a transit of HD 40307b can be considered as very 
unlikely, even if the star is actually significantly larger than de- 
duced from its parallax and bolometric luminosity. 

Figure 6 shows the IRAC photometry corrected for the 
baselines of the three time series and period-folded on the best- 
fitting Tq. The rms of this light curve is 143 ppm for a mean in- 
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Deduced parameters 
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Orbital semi-major axis a 
Orbital eccentricity e 
Argument of periastron a> 
M p sin i 

Orbital inclination i 
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Table 1. Parameters deduced for the three planets orbiting around HD 40307 from our MCMC analysis. Jump parameters are the 
model parameters that are randomly perturbed at each step of the MCMC. 




Fig. 5. Posterior probability density: time of minimum light To 
vs transit impact parameter b tr . 

terval of 1 .9 minutes, while the rms goes down to 29 ppm when 
the light curve is binned with an interval of 30 minutes. This il- 
lustrates well the high photometric potential that the cryogenic 
Spitzer had at 8 fim. Now that its cryogen is depleted, only the 
3.6 and 4.5 jjm channels remain operational, and similar pho- 
tometric precisions are not guaranteed, mainly because of the 
inhomogeneous intra-pixel sensitivity of the InSb detectors of 
IRAC (e.g. Knutson et al. 2008 and references therein). 



Fig. 6. Top: Light curve obtained after correction for the sys- 
tematics of the three IRAC light curves used in this work and 
folding on the best-fitting ephemeris. Bottom: same but after 
binning with an interval of 30 minutes. 

If present, the transit of a planet as small as 1.1 R @ (220 
ppm) should have been easily detected. To check this, we in- 
jected a fake transit of ~1 hour and 220 ppm in the Spitzer 
photometry, using the values reported in Table 1 for To and P. 
Our first MCMC analysis of these data was composed of two 
chains of 10 6 steps. Except in the 'burn-in' phases of the chains, 
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Fig. 7. Period-folded photometry with the best-fitting transit 
model obtained after insertion of a fake transit with a depth 
of 220 ppm and a duration of 1 hour. 



a transit was present in all the steps of both chains, indicating 
a firm detection. To estimate a formal transit probability, we 
performed two new chains for which no transit was allowed. 
For both models (with and without transit), the marginal like- 
lihood was computed from the MCMC samples like described 
by Chib & Jeliaskov (2001). The deduced Bayes ratio (e.g., 
Carlin & Louis 2008) was 1.4 x 10 23 in favor of the model 
with a transit, confirming the very firm transit detection. Fig. 
7 shows the obtained best-fitting transit model superimposed 
on the period-folded IRAC photometry corrected for the sys- 
tematics. The deduced transit depth is 232 + 23 ppm, in good 
agreement with the actual value of 220 ppm. 

We investigated the origin of the 'transit' marginally de- 
tected in the first IRAC time series alone. It appears that it was 
most probably due to the baseline model used in our prelimi- 
nary analysis. To recall, DD fitted a line in the light curve after 
cutting off the first three hours, while MG cut off the first hour 
and assumed a quadratic logarithmic function of time for the 
ramp model, i.e. the same function than shown in Eq. 1 but 
without the a\dt term. Figure 8 shows the photometry obtained 
after correction by the best-fitting baseline model obtained by 
least-square minimization, using Eq. 1 for the model, with and 
without the a\dt term. The residual correlated noise visible in 
the time series obtained when the a^dt term is not used could be 
easily 'compensated' by a shallow transit, while the correlated 
noise is nearly negligible when the a\dt term is used. Using 
only the first IRAC light curve and the HARPS measurements 
as data, we performed two new MCMC analysis, one with and 
one without the a\dt term in the baseline equation. The stellar 
radius was set to 0.86 R e . Without the a\dt term, the deduced 
transit probability was 97.1%, while it dropped to 0.9% with 
it. Comparing both models, the resulting Bayes ratio was 7.4 
in favor of the model with the a\dt term. We conclude thus 
that the significance of our transit detection in the first IRAC 
light curve was in fact smaller than estimated in our prelimi- 
nary analysis. This illustrates well how careful one should be 
when searching for very shallow transits in photometric time 
series presenting systematic effects of much larger amplitude, 
even if these systematics are rather well characterized. 
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Fig. 8. Light curves for the first run corrected for the baseline 
model and binned with a five minutes interval, without (top) 
and with (bottom) the linear term a\dt in Eq. 1. 

5. Conclusions 

We have used Spitzer and its IRAC camera to search for the 
transit of the super-Earth HD 40307b. The transiting nature of 
the planet could not be firmly discarded from our first pho- 
tometric monitoring of a transit window because of the un- 
certainty coming from the modeling of the photometric base- 
line. To obtain a firm result, two more transit windows were 
observed and a global Bayesian analysis of the three IRAC 
time series and the HARPS radial velocities was performed. 
Unfortunately, any transit of the planet during the observed 
phase window is firmly discarded, while the probability that 
the planet transits but that the eclipse was missed by our obser- 
vations is nearly negligible (0.26%). 
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